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Abstract — Protocol sequences are used for channel access in 
the collision channel without feedback. Each user accesses the 
channel according to a deterministic zero-one pattern, called the 
protocol sequence. In order to minimize fluctuation of throughput 
due to delay offsets, we want to construct protocol sequences 
whose pairwise Hamming cross-correlation is as close to a 
constant as possible. In this paper, we present a construction 
of protocol sequences which is based on the bijective mapping 
between one-dimensional sequence and two-dimensional array by 
the Chinese Remainder Theorem (CRT). In the application to the 
collision channel without feedback, a worst-case lower bound on 
system throughput is derived. 

Tags: Protocol sequences, collision channel without feed- 
back, cyclically permutable constant-weight codes, optical 
orthogonal codes. 

I. Introduction 
A. Background and Motivation 

OQ Randomness is commonly used in the design of multiple- 
access schemes. For example, in slotted ALOHA, each user 
transmits a packet with probability p independently. Imple- 
mentations of such random access schemes in practice usu- 
ally substitute random variables by pseudo-random numbers. 
However, high-quality pseudo-random number generation may 
be too complicated for applications, such as wireless sensor 
networks, where computing power is limited. The objective 
of this paper is to construct binary pseudo-random sequences, 
called protocol sequences, that are tailored to the quality-of- 
service requirements of interest, such as high throughput and 
bounded delay. 

Protocol sequences are used in multiple-access in the colli- 
sion channel without feedback (TJ. In this paper, we consider 
a time-slotted system, consisting of a number of transmitters 
and one receiver. A user sends a packet within the boundaries 
of a time slot. If exactly one user transmits in a time slot, 
the received packet is received successfully. If two or more 
users transmit in the same time slot, a collision is incurred, 
and the received packet is assumed unrecoverable. If no user 
transmits in a time slot, that time slot is idle. Since there is no 
feedback from the receiver, collision resolution algorithm such 
as the stack algorithm is not possible. Each user repeats his 
assigned binary protocol sequence periodically, and transmits a 
packet if and only if the value of the protocol sequence at that 
time slot equals one. We note that the transmission schedule is 
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independent of the data being sent and there is no cooperation 
among the users. 

The capacity of the collision channel without feedback is 
characterized by Massey and Mathys in (TJ. It is shown that the 
zero-error sum-capacity is e _1 . They use protocol sequences 
having the special property that the Hamming cross-correlation 
is independent of relative delay offsets. In fact, the Hamming 
cross-correlation they considered is a generalized notion of 
Hamming cross-correlation, which is defined for all nonempty 
subsets of users, not just for pairs of users. Protocol sequences 
with this property are called shift-invariant sequences J2). 
Shift-invariant protocol sequences have the advantage that 
there is no fluctuation in throughput no matter what the delay 
offsets are, and hence have the largest worst-case system 
throughput. Constructions of shift-invariant protocol sequences 
are considered in JT]— iEl- Nevertheless, shift-invariant pro- 
tocol sequences have the drawback that the period grows 
exponentially in the number of users J2J. Even if we relax 
this requirement and consider protocol sequence sets with 
only pairwise Hamming cross-correlation being constant, it is 
shown in |4] that the period grows exponentially in the number 
of users as well. Long period length has the disadvantage 
that individual or system throughput is invariant only if it 
is averaged over a long period. Individual users may suffer 
short-time starvation. In order to achieve short period length, 
we must seek for protocol sequences with some small variance 
in Hamming cross-correlation allowed. 

After the seminal work of Qj, more constructions of proto- 
col sequences are given in [5|-[9], sometime under the name 
of cyclically permutable constant-weight codes, or optical 
orthogonal codes. The main difference between these works 
in the literature and our construction is that, the protocol 
sequences in these papers are required to have small Hamming 
cross-correlation and auto-correlation. In our construction, 
Hamming auto-correlation may be very large. 

Another class of protocol sequences, called wobbling se- 
quences ifTOl . has period equal to M A , where M is the number 
of users, with worst-case system throughput provably larger 
than a positive constant that is approximately equal to 0.25 
when M is large. The result in this paper improves upon 
the wobbling sequences by constructing protocol sequences 
of order 0(M 2 ) or 0(M 3 ), depending on the models of 
user activity, while the guarantee of the worst-case system 
throughput remains the same. 

The construction proposed in this paper is based on the Chi- 
nese remainder theorem (CRT), and thereby the constructed 
sequences are called CRT sequences. We remark that the use 
of CRT in the construction of protocol sequences is not new. 
The constructions of optical orthogonal codes in [5|-[7| also 
employ CRT. However, the specific construction proposed here 
is novel. Although the analysis of CRT sequences is quite 
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involved, the generation of such sequences requires no more 
than computing linear operations in modular arithmetics 

The proposed CRT sequences have two special features, 
user-identification and frame- synchronization capability. The 
sender of each successfully received packet can be identified 
by looking at the channel activity only, without looking into 
the packet contents. Also, we can determine the start time of a 
protocol sequence, again, based on the channel activity only. 
One potential application based on this property is tracking 
targets by detecting energy pulses coming from multiple 
distributed sources, for example, as in a multistatic radar 
system or in a ultrasound based sensor network. While it may 
be difficult to identify the transmitting source of a single pulse, 
if each source employs a CRT sequence and sends multiple 
pulses according to the corresponding protocol sequence, then 
the identifiability property ensures that a detector can identify 
the transmitting sources of the pulses and make use of the 
information for more accurate tracking of the targets. 

B. Main Results 

The main construction is given in Section [XT] One of 
the key ingredients in the construction is the one-to-one 
correspondence between one-dimensional sequence and two- 
dimensional arrays via CRT. The correlation property of the 
constructed sequences is analyzed in Section Hill 

Applications to the collision channel without feedback are 
addressed in Section [TV] with two different user activity 
models ifTTIl . In the first model, an infinite backlog of data is 
assumed and the users are active throughout the transmission. 
We show in Theorem [7] that under this user activity model, 
we can achieve a throughput of 0.25 with 0(M 2+e ) sequence 
period, where M is the number of users. 

In the second model, users become active only if they have 
data to send, and remain idle otherwise. If a user becomes 
active, it is required that the user remains active for at least 
one period of the protocol sequence assigned. We show in the 
second part of Section HVl that if the number of active users is 
no more than one half of the number of potential users, and if 
the protocol sequences are sufficiently long, namely 0(Ajf 3 ), 
the receiver can detect the set of active users and determine 
their starting time correctly, even without any packet header. 

To facilitate comparison between different protocol se- 
quences, we introduce in Section [V] a notion called e- 
uniformity, which measures the variation of Hamming cross- 
correlation. Other applications of CRT are given in Section [VTI 

C. Definitions and Notations 

Let Z n be the ring of residues mod n for a positive integer n. 
We will reserve the letter L for sequence length. 

Definition 1. The components in a sequence of length L are 
indexed from to L — l. The time indices {0, 1, — 1} is 
identified with Z^. The Hamming weight of a binary sequence 
a(t) of length L is denoted by w a - The duty factor (T) of a{t) 
is the Hamming weight divided by the length, 

1 
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For two binary sequences a(t) and b(t) of length L, their 
Hamming correlation function is defined as 

L-l 

H ab (T) := a{t)b{t - r), 
t=o 

where r is the delay offset, and t — r is the difference in 
modulo-i arithmetic. When a(t) = b(t), H aa {r) is called 
the Hamming auto-correlation of a(t). When a(t) and b(t) 
are two different sequences, H a t,(r) is called the Hamming 
cross-correlation of a(t) and b{t). 

When the number of ones in a zero-one sequence is small 
in compare to the length, the sequence can be compactly 
represented by specifying the locations of ones. 

Definition 2. Given a sequence a(t) of length L, let the 

characteristic set of a(t), denoted by X a , be the subset of 
Zi such that t £ I a if and only if a(t) = 1. 
Shifting a sequence cyclically by r is equivalent to translating 
its characteristic set by r, with addition performed modulo L. 
Given a subset I in Z^, and r G Z^, we denote the translation 
of X by r as 

1 + t := {x + t e Z L : x <El}. 

Expressed in terms of the characteristic set, the Hamming 
correlation of sequences a(t) and b(t) equals 

H ab (T) = |X a n(I 6 + r)| 

for all t = 0, 1, . . . , L — 1, where |5| denote the size of a 
set S. 

II. The CRT Construction 

A. The CRT correspondence 

We shall construct sequences with length L — pq, where p 
and q are relatively prime integers. In subsequent discussions, 
we take p to be a prime number and q an integer not divisible 
by p. 

Define a mapping from Z p9 to the direct sum 

Gp t q '■= Zp © Zg 

by 

<& Pi q(x) := (x mod p, x mod q). 

By Chinese remainder theorem lfl2l p. 34], g is a bijective 
map. We will call $ p g the CRT correspondence. When the 
values of p and q are clear from the context, we write 
instead of $ p , q (x). 

It can be easily checked that the CRT correspondence is a 
linear map, meaning that 

<t>p,q( X + X ') = ®P,q( X ) + ®P,q( X ')- 

Here, the addition on the left hand side is the addition in Z pg , 
and the addition on the right hand side is the addition in G Ptq . 

A sequence s(t) of length L is associated with apxq array 
S(ii , t 2 ), where ti and t 2 range from to p — 1 and to q — 1 
respectively, via the relation 

S(t mod p, t mod q) = s(t). 
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The corresponding characteristic set of s(t) in Zl is mapped 
to a subset of G P:q via $ p g as well. 

Under the CRT correspondence, a one-dimensional cyclic 
shift of s(t) by one time unit is equivalent to a column- 
wise shift followed by a row-wise shift. One-dimensional 
correlation properties can be translated to the two-dimensional 
ones. 

For illustration, consider the time indices (0, 1,2,..., 14) 
as an integer sequence of length L = 15, p = 3 and q = 5. 
By the bijection $3,5, this integer sequence is mapped to 





10 

5 



G 
1 

11 



12 
7 
2 



3 

13 



9 
4 
14 



When (0, 1, ... , 14) is cyclically shifted to 
(14, 0, 1, 2, ... , 13), the corresponding array is 

"14 5 11 2 8' 
9 6 12 3 
4 10 1 7 13 

Note that the second array can be obtained from the first one 
by cyclically shift downward by one row and then to the right 
one column. 

For r = (ti,T2) £ Gp, q , define the Hamming correlation 
between two 2-dimensional arrays A and B by 

H ab (t) :=£)A(t)B(t-T), 
t 

with the subtraction calculated in G p<q , and t running over 
all elements in G p<q . It is easy to check that this definition 
of Hamming correlation is compatible with the 1 -dimensional 
analog, i.e., 

#ab(ti,T2) = Hab(r) 

with T\ = t mod p and Ti = r mod q. 

The Hamming correlation between two 2-dimensional ar- 
rays A and B can be expressed in terms of their characteristic 
set Za and I-q as 

^ab(t) = |i A n(x B +T)|, 

where 2a denotes the characteristic set of A, 

I A :={{i,j)£Z p ®Z q : A(i,j) = l}, 
and the addition in 2b + t is performed in Z p © 1 q . 

B. The CRT sequences 

We will construct sequences by specifying characteristic sets 
in G p _g. Rows and columns of matrices and arrays will be 
indexed by {0, 1, . . . ,p— 1} and {0, 1, . . . , q— 1} respectively. 

Definition 3. Let p be prime and q be an integer not divisible 
by p. For g £ Z p , we define 



{(g,l)t£G M :0<t<q}. 



(1) 



The notation (g, l)t simply means the sum of t copies of (g, 1) 
in G„ 



*p,q, 



(g,l)t:= (5,1) + (ff,l) + ... + (g,l). 



L 0,3,5 



-1,3,5 : 
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2,3,5 • 
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Fig. 1. The three arrays associated with characteristic sets Xo,3,5, Zl,3,5> 
3-2,3,5- 



We say that X g . p . q is generated by g, and g is the generator 
of I g ,p,q. If p and q are clear from the context, we write X g 



instead of 2, 



g,p,g- 



The elements in 1 g , p , q form an arithmetic progression in 
G Pi q with common difference (g, 1). We can re-write X g>P , q 
in the following form 

Zg.P.g = {(ffM) 6 G Pi q : < t < q}, 

with the product gt in the first component reduced mod p. 

Remarks: For g = 0, 1, . . . ,p — 1, the array with l g , P , q as 
the characteristic set contains exactly one "1" in each column. 
For g ^ 0, each block of p consecutive columns form a 
permutation matrix. (Recall that a permutation matrix is a 
square zero-one matrix with exactly one "1" in each row and 
each column.) 

Definition 4. (CRT sequences) For g = 0, 1, . . . ,p — 1, define 

the CRT sequence generated by g, denoted by s g , p , 9 (i), be the 
binary sequence of length L obtained by setting 



1 if %, q (t) £ Xg,p, q 

otherwise. 



s g ,p,q(t) 

We will write s g (t) if the values of p and q are understood. 

3 and q = 5. The three characteristic sets 



Example 1. p 

are: 



2o,3,s = {(0, 0), (0, 1), (0, 2), (0, 3), (0, 4)}, 
Ti, 3 ,5 = {(0,0),(l,l),(2,2),(0,3),(l,4)}, 
22,3,5 = {(0,0), (2,1), (1,2), (0,3), (2, 4)}. 

The three arrays are shown in Fig. Q] The top left corner in 
each array is the (0, 0) -entry. The generated CRT sequences 
are listed as follows 



s (t) 

*2(*) 



10010 0100100100 
1111100000 00000 
10010 00100 01001. 



CRT sequences satisfy the following properties. 
Theorem 1. 

1) The Hamming weight of each CRT sequence is q. 
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2) so(t-\-p) — so(t), i.e., the least period of sequence So(t) 
is p. 

3) For each residue i mod q, there is exactly one "1 " at 
t = i,i + q,i + 2q, . . . ,i + (p — l)q in s g (t). 

Proof: As the characteristic set I g , P , q contains q elements 
and the CRT correspondence $ preserves Hamming weight, 
the first statement in the proposition follows immediately. 

For 3 = 0, the elements in characteristic set Io,p,q have the 
first coordinate identically equal to zero. For j = 0, 1, . . . , q — 
1, The pre-image of (0, j) under $ is a multiple of p. We have 
So(t) = 1 whenever t is a multiple of p. 

For the last statement in the proposition, consider the image 
of i + kq, for k = 0, 1, . . . ,p — 1, under the mapping $, 

$(i + kq) — (i + kq mod p, i mod q). 

The second coordinate is constant. If k goes through to 
p — 1, $(z + kq) will go through the ith column in the array 
with characteristic set I g:P: q. Because there is exactly one "1" 
in each column, there is exactly one "1" in positions t = 
i,i + q,...,i + (p-l)q. ■ 

III. Correlation Properties of CRT Sequences 

We continue to use the notation that p is a prime number, 
q is an integer not divisible by p. Let L — pq denotes the 
sequence length. In this section, we determine the Hamming 
correlation of the CRT sequences. In this section, we will use 
the notation x for the remainder of x after division by p. We 
distinguish the two additions mod p and mod q by (B p and 
® q , respectively. 

For g = 0, 1, . . .,p — 1, the p x q array corresponding to 
sequence s g (t) is denoted by A g . Recall that the characteristic 
set of A g is given in (Q3, and A g (i,j) = 1 if and only if 
i = jg mod p, for i = 0, . . . ,p — 1 and j = 0, . . . , q — 1. For 
notational convenience, we let 

H gh (Ti,T 2 ) := H AfjAh (n , r 2 ) . (2) 

A. Hamming Cross-correlation 

Let g and h be two distinct elements in Z p . As argued in 
the previous section, the cross-correlation of CRT sequences 
s g (t) and s/j(i) is equivalent to the cross-correlation of the 
associated p x q arrays A g and A/j, namely by counting the 
number of elements in common to I g and X/j + (t\, t 2 ). By 
definition, the translation Z/, + (n, r 2 ) of I/j is 

{((jh)® p n,j® q T 2 ) : j = 0,1,..., (3) 
By a change of variable, (O can be written as 

{(((j e q r 2 )h) (B p Ti,j) : j = 0,1,..., q-1}. 

where Q q denotes subtraction mod q. In the calculation of 
(j Qq T 2 )h, j& q T2 must be first reduced to an integer between 
and q — 1 before we reduce modulo p. The intersection of 
I g and Xh + (ti,T2) equals to the number of solutions to 

xg = [[xQ q T 2 )h) (Bp ti mod p. (4) 

for x = 0, 1, . . . , q — 1. 



For nonzero h, we can divide both sides of © by h and 
re-write it as 

x{h~ 1 g) = (x Qq t 2 ) (B p (/i _1 ri) mod p. 

For each fixed t 2 , as t\ runs through Z p , hT X T\ also runs 
through Z p . Therefore, the distribution of Hamming cross- 
correlation between s g (t) and Sh(t) is the same as the 
distribution of Hamming cross-correlation between s g / h (t) 
and s\(t). From this observation, it suffices to take h = 1 
without any loss of generality. We will consider the Hamming 
cross-correlation between s g (t) and si(t), for g — and 
.9 = 2,3, ...,p- 1. 

The following simple lemma is used repeatedly in the 
derivation of Hamming cross-correlation properties. 

Lemma 1. For each b £ Z p , the number of solutions to 

x = b mod p 

for x going through d consecutive integers c, c+1, . . . c+d—1, 
equals 

{d/p if p divides d, 

[d/p\ + 5 otherwise, 

where S equals either or 1. 

Proof: Suppose that d is divisible by p. If we reduce the 
integers c, c + 1, . . . ,c + d — 1 mod p, we have each element 
in Z p repeated d/p times. Hence, for each b £ Z p , there are 
exactly d/p integers in {c, c+ 1, . . . , c+d— 1} whose residue 
mod p equal b 

Now suppose that d is not divisible by p, we divide the d 
consecutive integers into two parts. Among the first [d/p\p 
integers, for each b £ Z p , exactly [d/p\ equals b mod p. The 
residues of the remaining d — [d/p\p integers are distinct. The 
number of integers in {c, c + 1, . . . , c + d — 1} whose residues 
equal b is either [d/p\ or [d/p\ + 1. ■ 

We obtain the following theorem immediately from 
Lemma Q] 

Theorem 2. The Hamming cross-correlation of s\(t) and 
So(t) is equal to either \q/p\ or \q/p\ + 1. 

Proof: If we put g — 1 and h = in (0|, we get 

x = Ti mod p. 

The number of integers in {0, 1, . . . , q— 1} that equal t\ mod p 
is either [q/p\ or |_<?/pj + 1 by Lemma Q] ■ 
From now on, we assume q > p, which is the case of 
practical interest. 

Theorem 3. Suppose q > p. Let m be the quotient of q divided 
by p, i.e., m = \q/p\, and let g £ Z p , ^ g 1. Let q be 
the residue of q mod p, and 

b g = (g - iy X q mod p. (5) 

The Hamming cross-correlation between s g (t) and Si(t) is 
bounded between 

rn — 1 and to + 1 if < b g < p — q, or 
m and m + 2 if ' p — q < b g < p. 
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The proof of Theorem [3] is in Appendix lAl We remark that 
in Theorem [3] b g is neither nor p — q. The value of b g is 
nonzero mod p because both g — 1 and g are nonzero. On 
the other hand, b g is equal to p — q only if g = 0, which is 
excluded by assumption. 

With more careful book-keeping, we can determine exactly 
the frequencies of occurrence of Hamming cross-correlation. 

Definition 5. Let the distribution of H g x(rx,T2) be 

N g (j) := |{(ri,r 2 ) g G pq : H gl ( n ,T 2 ) = j}\, (6) 

for j = 0,1,...,?. 

Theorem 4. Wi'f/z notation as in Theorem\3\ we have 

1) N {m) = (p- q)q, N Q (m + 1) = gg. 

2) If < b g < p - q, then 



N g {m — 1) = rj 

N g (m) = g(p -q)-2r) 
N g (m + 1) = qq + r), 



(7) 
(8) 
(9) 



where 



i] := mb g (p - b g 



3) If p — q < bg < p, then 



N g (m) = q(p-q) + > 
N g {m + 1) = qq-26, 
N g (m + 2) = 9 



(10) 

(11) 
(12) 



where 



9:= (m + l)(p-b g )(q-\ 
The proof is relegated to Appendix 151 



•P). 



Example 2. For the CRT sequences in Example [T] with 
parameters p = 3 and q = 5, we tabulate the distribution 
of Hamming cross-correlation between s s (t), and si(i), for 
g 7^ 1, as follows. 
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^(2) 


^(3) 
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5 


10 
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2 


7 


6 


2 



We note that if g = 2, then b g is 2 mod 3. By Theorem [3] 
H2i(t) equals 1, 2, or 3. The Hamming cross-correlation is 
distributed according to the third part of Theorem |4] 

When q = ±1 mod p, the Hamming cross-correlations have 
three distinct values, for all pairs of distinct CRT sequences 
chosen from si(t), S2{t), . ■ . , s p _i(i). 

Theorem 5. Suppose q > p. 

1) Let q be of the form mp + 1. For g = 2, 3, . . . ,p — \, 
H g i (r) is between m — 1 and m + 1. 

2) Let q be of the form mp+(p— 1). -For <? = 2, 3, . . . ,p— 1, 

(r) k between m and m + 2. 



Proof: For the first part of the theorem, we have q equal 
to 1 mod p. So 

= (.9 - I)" 1 ? = (5 - I)" 1 mod p. 

When g runs from 2 to p — 1, the value of g — 1 runs over all 
elements in Z p except and p— 1. For odd p, the multiplicative 
inverse of p — 1 is itself. Hence the range of (<? — mod p 
is Z p \ {0, -1}. We thus obtain < b g < p - 1. The result 
now follows from Theorem [3] 

The second part can be proved similarly. ■ 

B. Hamming Auto-correlation 

The distribution of the Hamming auto-correlation can also 
be determined explicitly. 

Theorem 6. 



H 00 (t) 



q if t is a multiple of p, 



otherwise. 
For g — 1, 2, . . . ,p — 1, and k = 0, 1, . . . , q — 1, 
'g-fc ; /$(r) = ±(.g,l)fc, 







otherwise. 



Proof: The first statement follows by the second part of 
Theorem Q] that s Q (t) has least period p. 

For the second statement, we note that (g, l)t runs through 
the whole g roup Gp^q us i runs from to pq — 1. This follows 
from CRT, because for any x and y, 

gt = x mod p 
t = y mod g 

has one and only one solution mod pg. Via a bijective 
mapping t H> (<?, the set I g can be regarded as an arith- 
metic progression {0, 1, 2, . . . , g — 1} in Z pg . The intersection 
I g n (I g + is a. subset of I g that is also an arithmetic 

progression. The intersection I g n (I g + is nonempty if 

and only ±(i, j) 6l s , If |I g n (X g + is nonzero, then 

is equal to ±(g, l)k for some = 0, 1, . . . , g — 1, and 

\I g n(I g + k(g,l))\=q-\k\. 



IV. Application to the Collision Channel without 
Feedback 

Consider a time-slotted collision channel with K trans- 
mitters and one receiver as described in the introductory 
section in this paper, and assume that there is no cooperation 
among the users, and no feedback from the receiver. We 
regard a packet as a Q-ary alphabet chosen from the alphabet 
set Q = {1, 2, . . . , Q}. The channel output at slot t equals 

if no user transmits at slot t, 
* if two or more users transmit at slot t, 
x if exactly one user transmits a packet with content x. 
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Each user is assigned a deterministic and periodic zero-one 
sequence, called protocol sequence (fl. For i = 1, 2, . . . , K, 
the protocol sequence associated with user i is denoted by 
Sj(t), which is a periodic binary sequence of period L. As 
there is no feedback from the receiver and no cooperation 
among the users, each user has a relative delay offset r, which 
is random but remains fixed throughout the communication 
session. User i sends a packet at slot t if Si(t + r) = 1, and 
remains silent if Si(t + r) = 0. 

We consider two different factors in the system model: 
(i) systems with packet header or without packet header, 
and (ii) systems with permanent user activity or partial user 
activity. For system without packet header, the receiver has 
to identify the set of active users (user identification), and 
determine the sender of each successfully received packet 
(packet identification). If packet header is present, information 
such as user identity is obtained readily by looking up the 
header, and there is no user and packet identification problem. 

For system with partial user activity, the number of active 
users M is smaller than the total number of potential users T. 
When a user changes from inactive to active, packets are sent 
according to the protocol sequence assigned. It is assumed that 
after all the data are sent, the user must remain inactive for at 
least L slots before becoming active again. On the other hand, 
if permanent user activity is assumed, the number of active 
users M is equal to the total number of users, and each user 
sends packets periodically according to the protocol sequence. 



A. Permanent User Activity with Header 

For systems with a packet header, the users are simply 
identified by packet headers. We calculate the achievable 
throughput when CRT sequences are used. Let the number of 
users be M. We pick M sequences from the CRT sequences 
generated with parameters p and q. To take advantage of the 
three-valued property mentioned in Theorem [5] we pick a q 
which is —1 mod p. 

Construction 1. Let p be a prime number and k a positive 
integer. We choose q = kp — 1. The CRT construction in 
Definition [4] yields p protocol sequences of length L — kp 2 ~p 
and weight w = kp — 1. We pick any M sequences from this 
set of CRT sequences. 

Since the Hamming weight is kp — 1 and the pairwise 
Hamming cross-correlation is at most k + 1 by Theorem [5] 
the total number of successful packets (summed over all M 
users) is lower bounded by 



M[kp - 1 - (M - l)(k + 1)]. 
By completing square, we can write the above as 
k(p+l)\ 2 , / fc(p + l) 



(fc + 1) 



- M - 



2(fc + i)y V2(fc+i) 



(13) 



(14) 



Consider (fl4] > as a function of M. The maximum value is 
obtained when 

k(p + 1) 



and the maximal value is 

(p+1) 2 



(16) 



4 fc+1 

After dividing the above by the period, we obtain the following 
lower bound on throughput. 



Theorem 7. Let 



M = 



k(p+ 1) 



-2(fc + l) 

in Construction^ The system throughput is lower bounded by 

1 (p + 1) 2 k 2 fc + 1 
4 kp 2 — p fc+1 kp 2 — p 



(17) 



Proof: Consider the expression in ( fT4b as a function of M, 
and denote it by f(M). The coefficient of M 2 is —(k+ 1). If 
we take M = M, we deviate from the optimal value of M* 
by at most 1. Therefore, 

f(M*) - f(M) = (fc + 1)(M - M*) 2 < k + 1. 

If we put M = L 2(fc+ 1 j -I m Q~U' we ^ ro P f rom me optimal 
value in ( fTBI l by at most fc + 1. Hence the total number of 
successful packets, divided by the period, is lower bounded 
by CCD. ■ 
Theorem [7] provides a lower bound on the worst-case 
throughput. The mean system throughput, averaged over all 
delay offsets, is however much higher than the lower bound. 

Example 3. We consider an example with M = 19 users, 
using CRT sequences with p = 37. The throughput is plotted 
against sequence length with increasing fc, while keeping the 
duty factor fixed at 1 /p. We compare the lower bound in (TTTb 
of worst-case throughput with the average throughput obtained 
by simulation in Fig. [2] For each fc, 100000 delay offset combi- 
nations are randomly generated. Beside the mean throughput, 
the maximum and minimum throughput obtained among these 
100000 delay offset combinations are also plotted. We can 
observe that the mean system throughput is about 0.31, in 
accordance with the theoretical value 
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= 0.314. 
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When the length increases, the lower bound approaches 
0.25(p+ l) 2 /p 2 - l/p 2 = 0.263. 

Asymptotically, if we increase fc and p in such a way that 
fc increases much slower than p, we obtain a lower bound on 
the system throughput of 1 /4. 

Theorem 8. For arbitrarily small t,8 > 0, there exists an 
infinite class of protocol sequence sets of length 0{M 2+<L ), 
where M is the number of sequences, with system throughput 
lower bounded by 0.25 — 5 for all sufficiently large M. 

Proof: We choose fc = log(p). The lower bound of system 
throughput in Theorem [7] tends to 1 /4, with 

L~kp 2 - fc(2A/) 2 - AM 2 logM. 

We can find an integer Mo sufficiently large such that 
4Mq logM is less than Mp +e . The theorem then holds for 
all M > M . M 
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Sequence Peirod 

Fig. 2. System Throughput of CRT Sequences for 19 Users, p = 37. 

B. Partial User Activity without Header 

In this application, we use a modified version of the CRT 
correspondence. Let 7 be the multiplicative inverse of p in Z ? , 
i.e., 7p = 1 mod q. Since p and g are relatively prime, such 
inverse exists. Define the mapping &' p : Z^^ — >• G p g by 

^p.q( x ) '■— ( x m °d p, 72; mod g). 

It can also be shown that &' p q (x) is a group isomorphism 
between Z pq and G P:q . 

Because of its unfavorable Hamming auto-correlation prop- 
erty, the sequence generated by g = is not used in this 
application. Given a prime number p, the system supports p—1 
potential users. We label the users from 1 to p — 1. 

Construction 2. Let p be a prime number, and g be an integer 
relatively prime to p. Choose 7 as described above. For g = 
1, . . . ,p — 1, construct the CRT sequence s g (i) of length pq 
by setting 




1 if <!>;„/:/:, /,, 
otherwise. 



The sequence generated by g is assigned to user g. We have 
used <£>p q (x) in the above construction, instead of <& p , q (x) as 
in Construction Q] The cross-correlation properties of the re- 
sulting modified CRT sequences are exactly the same as in the 
previous section, because the proof in the previous section is 
essentially about two-dimensional Hamming cross-correlation. 
If $ is replaced by all theorems about Hamming cross- 
correlation and auto-correlation also hold. 

Sequence synchronization 

For systems without a header, we need to identify the sender 
of a successful packet. Also, as users may come and go, we 
also need to determine when a user becomes active. We show 
that these two tasks can be achieved by merely observing the 
channel activity, that is whether a time slot is idle, containing 



a collision or a successful transmission, without looking into 
the packet contents. 

Definition 6. For each time index t, let c(t) be if it is an 
idle slot, 1 if exactly one user transmits a packet, or * if 
two or more users transmit. We call c(t) the channel-activity 
signal. We say that c(t) is matched to Si(t) at time to if 
Vt = 0, 1, . . . , L - 1, Si(t) = 1 => c(t + 1) = 1 or *. 

The receiver stores the channel-activity signal in a first-in- 
first-out queue. 

We want to determine (a) the time when a user becomes 
active, and (b) the time when an active user change status 
from active to idle. The receiver keeps track of the active 
users by maintaining p — 1 Boolean variables active(i), for 
i = 1, 2, . . . ,p — 1. The value of active(i) is set to FALSE if 
the user is idle, and TRUE if the user is active. To determine 
whether user i becomes active at time to, the receiver wait 
until time to+L. At that time the channel-activity signal c(to), 
c(to + 1), . . . , c(to + L — 1) are available. The receiver then 
checks whether c(t) is matched to Si(t) at time to- If there is 
a match, we declare that user i is active and the starting time 
of user i is stored in variable start(i). We summarizes the 
procedure in Algorithm [T] below. 

Algorithm 1 Determining when a user becomes active or 
inactive at time to. 

l: for i = 1, 2, . . . ,p — 1 do 

2: if active(i) = FALSE then 

3: if c(t) is matched to Sj(i) at to then 

4: active(i) «- TRUE 

5: start(i) <— to 

6: end if 

7: else 

8: if to — start{i) is a multiple of L and c(t) is not 
matched to Sj(t) at time t then 

9: active(i) «- FALSE 

10: start(i) <- 

li: end if 
12: end if 
13: end for 



Under some conditions on the number of active users and 
period length, we can show that the above algorithm is able 
to identify the starting time of each user. 

Theorem 9. Let p be a prime number, and q > 2p 2 be an 
integer relatively prime to p. Construct p—\ CRT sequences by 
Construction [2] and assign them to users 1 to p — 1. Suppose 
that at most (p + l)/2 users are active at the same time. 
Then Algorithm Q] can successfully identify the active users 
and determine their starting time. 

The proof of Theorem [9] is given in Appendix O 

Example 4. Consider the CRT sequences generated with 
parameters p = 7, q = 8 and 7 = 7. The period is equal to 
56. Suppose that users 1, 2, 3, 4 and 6 are active. The relative 
delay offsets of users 1 and 2 are 0, and the relative delay 
offsets of users 3, 4 and 6 are 1. The CRT sequences si(t) 



x 



si(t) : 10001000 00000001 00010000 00000010 00100000 00000100 01000000 10 

s 2 (t) : 10010000 00000100 00000001 00100000 00001000 00000010 01000000 10 

s 3 (t - 1) : 1000001 00000100 00010000 00000000 10000010 00001000 00100000 1 

s 4 (t - 1) : 1000010 00010000 00000001 00001000 00000000 10000100 00100000 1 

s 6 (t - 1) : 1000000 00000000 00000000 00000000 00000000 00000000 00111111 1 1 

c(t) : **011011 00010*01 000*000* 00101010 10101010 10001*10 0**11111 ** 

s 6 (t) : 10000000 00000000 00000000 00000000 00000000 00000000 01111111 

Fig. 3. CRT sequences in Example [4] Sequence periods are indicated by underbraces. 



to S4 (t) and s 6 (t) and the induced channel-activity signal c(t) 
are shown in Fig. [3] 

By comparing with Se(t), the receiver declares that the 
channel-activity signal is matched to Se(t) at time 0. The 
receiver cannot distinguish whether user 6 starts transmitting at 
time or 1, and erroneously detects that user 6 becomes active 
at time 0, while the actual delay offsets is 1. Thus AlgorithmQ] 
fails in this case. 

The preceding example illustrates that Algorithm Q] does not 
work if there are too many active users. Theorem|9]asserts that 
such error in synchronization does not occur if q > 2p 2 and 
the number of simultaneously active users does not exceed 
(p + l)/2. 

If q = ±1 modp, the Hamming cross-correlation of the 
resulting CRT sequences is three-valued (see Theorem[5]l. For 
these preferred choices of q, we can improve Theorem [9] by 
relaxing the requirement q > 2p 2 to q > p 2 . 

Theorem |9f. Let p be a prime number, and q > p 2 be 
an integer such that q = ±1 mod p. Construct p — 1 CRT 
sequences by Construction [2] and assign them to users 1 to 
p — 1. Then Algorithm [7] can identify the active users and 
determine their starting time, provided that the number of 
simultaneously active users is no more than [p + l)/2. 

The proof of Theorem |9| is analogous to the proof of 
Theorem [9] in Appendix ICl 

Erasure Correction and Throughput 

We have shown in Theorem [9] that the receiver is able to 
figure out the delay offsets of the active users. From a specific 
user's point of view, the channel reduces to an erasure channel. 

In the remainder of this section, we pick q to be an integer 
of the form kp + 1, for some integer k > p, so that q is larger 
than p 2 and Theorem |9j can be applied. 

For each user, the number of successfully received packets 
in a length is lower bounded by 

(fep+l)-(£±i-l)(fc + l). (18) 

The first term kp + 1 is the total number of packets sent by a 
user in a period. The factor k + 1 is the maximum Hamming 



cross-correlation from Theorem [5] After some simplifications, 
( [18] ) can be written as 

D(p,k):=^[kp + k-p + 3}. (19) 

An erasure-correcting code can be applied across the pack- 
ets in a period. We pick Q to be a power of prime such that 
Q > kp + 1. We encode D(p, k) information packets using 
a shortened Reed-Solomon (RS) code of length kp + 1 over 
the finite field of size Q. The kp + 1 encoded packets are 
sent out according to the assigned protocol sequence. Because 
of the maximal-distance separable property of RS codes, we 
can recover the information packets. We have thus proved the 
following 

Theorem 10. Suppose that there are (p + l)/2 active users 
out of p — 1 potential users in the collision channel without 
feedback, where p is an odd prime. Each active user can send 
D(p, k) information packets in a period of p(kp+l) time slots, 
where D(p, k) is given in (|19t and k is an integer larger than 
or equal to p. In particular, if we take k = p, the resultant 
CRT sequences have period p 3 +p, and when (p+l)/2 users 
are active, the system throughput is lower bounded by 

£+1 . o.sy+3) > 

2 p 3 + p ~ 



Example 5. We pick p = 19, k = 19 and q = p 2 + 1 = 
362 in Theorem [TOl Generate 18 CRT sequences of length 
pq = 6878 by Construction Q] Pick Q = 512 = 2 9 , which is 
a prime power larger than the Hamming weight of the CRT 
sequences under consideration. Using a shortened RS code of 
length 362 and dimension £>(19, 19) = 182 over the finite 
field of size 512, we encode 182 information packets in each 
period for each active user. When 10 users are active, the 
total number of information packets sent through the system 
is 182 x 10 = 1820, achieving system throughput no less than 
1820/6878 = 0.265. 
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V. Comparison with Other Protocol Sequences 

In order to compare the variation of Hamming cross- 
correlation due to delay offsets, we introduce in this section a 
measure of deviation called e-uniformity. 

Given any two binary and periodic sequences a(t) and 
b(t), let E r [i? b(r)] be the expectation of Hamming cross- 
correlation, with delay offset r chosen uniformly at random 
over a period. 

Definition 7. We say that the Hamming cross-correlation H ab 
is e-uniform if 

\H ab (v)-E T [H ab ( T )] 



(20) 



E T [H ab { T )] 

for allT = 0, 1,...,L — 1. A sequence set is called e-uniform 
if e is the smallest number such that each pair of distinct 
sequences is e-uniform. We say that a sequence set is pairwise 
shift-invariant if it is 0-uniform. 

In other words, H ab (r) is e-uniform if for all delay offsets 
t, the percentage difference between H ab (T) and the mean 
is between — e and e. The notion of e-uniformity is the 
same as the normalized distance between Hamming cross- 
correlation and the expected value. 

The sum of the Hamming cross-correlation over all relative 
delay offsets in a period equals w a w b , where w a and w b 
denote the Hamming weight of a(t) and b(t) respectively (cf. 
Lemma [3] in Appendix |B). If we take the average over all 
delay offsets r, then 

. », W n Wh 

E T [H ab ( T )\ = 

Hence, the definition of e-uniformity in ( f20b can be written as 

\H ab (v) - w a w b /L\ 
w a w b /L ~ 

A lower bound on the worst-case throughput similar to ( [T7l i 
can be expressed in terms of e-uniformity. Suppose that there 
are K active users and each of them is assigned a sequence 
from an sequence set of length L and Hamming weight w. 
Suppose that the sequence set is e-uniform. By the union 
bound, a user can successfully send at least 

w-(K -l)(l + e)(w 2 /L) 

packets in each period. Individual throughput is lower bounded 
by 

f-(K-l)(l + e)f 2 , 

where / is the duty factor w/L. It can be easily seen that for 
fixed duty factor / and number of users K, a smaller e yields 
a larger lower bound on individual throughput. 

Example 6. (Constant- weight cyclically permutable codes [5 1) 
In (5] p. 948] Example 5, a protocol sequence set of period 
156 and Hamming weight 12 is presented. The number of 
sequences is 169. It is shown that the Hamming cross- 
correlation is between and 3, and the mean Hamming cross- 
correlation equals (12 2 )/156 = 12/13. The maximal deviation 
from the mean is 3 — 12/13 = 27/13. This sequence set is 
thus (27/12)-uniform. 



Example 7. (Prime sequences lfT3"l ) Given a prime p, we can 
construct a sequence set with period p 2 , Hamming weight p 
and Hamming cross-correlation no more than 2. The mean 
Hamming cross-correlation is 1. For each pair of distinct 
prime sequences, the maximal \H ab (v) — E T [H ab (T)]\ over 
all possible delay offset v is equal to 1. The prime sequences 
are therefore 1 -uniform. 

Example 8. (Extended prime sequences [14]) By padding 
p — 1 zeros after every "1" in a prime sequence, we obtain a 
sequence set with period p(2p — 1), Hamming weight p and 
Hamming cross-correlation either or 1 . The mean Hamming 
cross-correlation is pj (2p — 1), which is roughly equal to one 
half. We can check that the extended prime sequences are 1- 
uniform. 

Example 9. (Wobbling sequences iflOl ) Based on prime se- 
quences, a class of (l/p)-uniform sequence sets is constructed 
in IflOl . The number of sequences is p and the sequence period 

4 

IS p . 

Example 10. (Shift-invariant and pairwise shift-invariant se- 
quences Q— [4]) All pairwise shift-invariant is by definition 
0-uniform. It is shown in |2| that the period grows exponen- 
tially in the number of sequences. For pairwise shift-invariant 
sequences, the period is also shown to grow exponentially with 
the number of users |4|. 

The examples above are presented in an order such that 
the e-uniformity is decreasing. We can see that the sequence 
period increases as we go down the list. 

Theorem 11. The sequences obtained by Const ruction\l\with 
period 0(kp 2 ) are (1/k)- uniform. 

Proof: The Hamming cross-correlation in Construction [T] 
is obtained by Theorem [5] The mean Hamming cross- 
correlation is 

w 2 (kp — l) 2 kp — 1 
L p(kp — 1) p 

The maximal difference between Hamming cross-correlation 
and the mean is 



k± 1 - 



kp — 1 



1, 

= - ±P-1 
P 

< -(p+1) 
P 

kp - 1 



P 



p+1 
kp — 1 



kp 



P 



2 <> 



We thus have an 0(l/fc)-uniform sequence set of period 
0{kp 2 ). ■ 

The trade-off between period length and e-uniformity is 
summarized in Table U 

To compare with the wobbling sequences, we can take k 
roughly equal to p. This yields CRT sequences with are (1/p)- 
uniform and period 0(p 3 ). The e-uniformity is roughly the 
same as wobbling sequences but the period is shorter than the 
period of the wobbling sequences. 
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Period 


e 


Prime sequences 


P 2 


1 


Extended prime sequences 


p(2p - 1) 


1 


CRT sequences 


0{kp A ) 


0(l/k) 


Wobbling sequences 


P 4 


l/p 


Shift-invariant sequences 


exponential in p 






TABLE I 

Tradeoff Between Period Length and e-uniformity for Various 
Sequence Sets for p Users 



VI. Other Applications 

A. Extension to Multiple Data Rates 

To support service with multiple data rates, we need pro- 
tocol sequences with different duty factors; a sequence with 
larger duty factor is assigned to a user with higher data rate 
requirement. The CRT construction can be extended to cope 
with multiple data rates. For the sequence s g , p ^(t) generated 
by g, we replace the characteristic set I g , P , q by 

I g , p , q U (lg tM + (1, 0)) U . . . U (Ig tM + (k, 0)) 

for some positive integer k. We note that the above is a union 
of disjoint sets. The resulting sequence has duty factor k/p. 
The measure of uniformity however does not change; if the 
original CRT sequence set is e-uniform, the extended sequence 
set is still e-uniform. 

B. Application to Multi-Channel Network 

In this network, users can send data to each other, and we 
have a fully connected system topology. The total bandwidth 
is divided into p subchannels and each subchannel is assigned 
to a user. Let the subchannel assigned to user i be denoted by 
subchannel i, for i = 1,2, ... ,p. Each user is also assigned 
a CRT sequence. The half-duplex model is assumed, so that 
each user cannot transmit and receive at the same time. Users 
always receive in their assigned subchannel. In a period of L 
slots, user i can pick one user, say user j, and send packets 
to user j via subcarrier j using user i's CRT sequence. In 
one sequence period, user i either: (i) receives packets in 
subcarrier i for the whole period, or (ii) sends packets to user j 
in L/p packets using subcarrier j and receives packets in the 
remaining L — L/p packets in subcarrier i. 

The worst-case scenario occurs when user i is sending 
packets to some other user, and all the remaining users want 
to send packets to user i in the same period. Using CRT 
sequences, we can show as in the multiple-access case that 
the worst-case throughput between each pair of users is lower 
bounded by a positive constant. 

VII. Conclusion 

A class of protocol sequences, called CRT sequences, whose 
Hamming cross-correlation is highly concentrated around the 
mean value is given. When CRT sequences are applied to 
the collision channel without feedback, we obtain a tradeoff 
between worst-case throughput and the sequence period. The 
generation of CRT sequences involves only simple modu- 
lar arithmetics, and provides a low-complexity solution to 
multiple-accessing in wireless sensor network. 



Appendix A 
Proof of Theorem[3] 

In this appendix, g denotes an element in Z p \ {0, 1}, and 

a g (n,T 2 ) := (g - 1) -1 (ti - f 2 ) mod p. (21) 

As defined in (0, b g denotes (g — l)^ 1 ^ mod p. We will also 
use the indicator function I defined as 

(l if Pis true, 
[0 if Pis false. 

Recall that x stands for the remainder of x after division by p. 

Lemma 2. The Hamming cross-correlation H g i(Ti,T 2 ) as 
defined in (0 satisfies the following properties: 

1) H g i{Ti,T2) equals the number of solutions to 

x = a g (Ti,T 2 ) ©p b g l(0 < x < r 2 ) mod p, (22) 

for x = 0, 1, . . . , q — 1. 

2) Let {t\,t 2 ) and (t[,t 2 ) denote two {2 -dimensional) 
relative delay offsets. If t\ — t[ and r 2 = r 2 + kp 
for some integer k, then 

H gl (n,r 2 ) = H gl (T{,4). 



Proof: After setting the value of h in to 1, we obtain 

xg = (x Q q r 2 ) ffi p t\ mod p. (23) 

We want to show that the number of solutions to (F23l . for 
x = 0, 1, . . . , q — 1, is the same as the number of solutions 
to d22l) . 

Consider x in two disjoint ranges: (i) < x < r 2 , and 
(ii) i~2 < x < q. In the first case, x Q q t 2 is congruent to 
x + q — t-2 mod q. So, for < x < r 2 , d23l is equivalent to 

xg = x (Bp q Q p f 2 ® p T\ mod p (24) 

where q and f 2 are residues of q and r 2 in Z p , respectively. 

In the second case, for x = t 2 , t 2 + 1, • • • , q — 1, (l23T > is 
equivalent to 

xg = x Op t 2 (Bp n mod p. (25) 
We combine d24b and d25b in one line as 

x(g — 1) = — r 2 © p rj © p ql(0 < x < r 2 ) mod p. 

Since g is not equal to 1 by assumption, we can divide by 
(g — 1) and obtain (f22j). This proves the first part of the lemma. 

The second part of the lemma is vacuous if q < p. So we 
assume q > p. (The case q — p is excluded because it is 
assumed that q is relatively prime with p.) It is sufficient to 
prove the statement for t% = t[ and t' 2 = r 2 + p, namely, the 
number of solutions to ( 122b and the number of solutions to 

x = a g (r 1 ,T 2 + p) + b g I{Q < x < t 2 + p) mod p (26) 

for x = 0, 1, . . . , q — 1, are the same. We note that a g (ri, t 2 ) 
is equal to a ff (ri,r 2 +p). However, the arguments inside the 
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indicator function are different. We divide the range of x into 
three disjoint parts: 

Xi := {0,1,- ..,75-1}, 

X 2 ■= {t 2 ,t 2 + 1, ...,r 2 + p- 1}, 

X3 ■= {r 2 +P,t 2 +p+ 1, .. .,q - 1}. 

For x G X\, 

1(0 < x < t 2 ) = 1(0 < x < t 2 +p) = 1. 

Therefore ( f22T > and d26*| i have the same number of solutions 
for x in X\. For a; G X 2 , both (l22l and (l26b have exactly one 
solution by Lemma Q] For x G X3, we have 

1(0 < X < T 2 ) = 1(0 < X < T 2 + p) = 0, 

and hence (l22l and ( 126*1 ) have the same number of solutions 
for x G X3. In conclusion, the number of solutions to ( f22l 
and ( f26l for x E XiU X 2 \J X3 are the same. This finishes the 
proof of the second part of the lemma. ■ 
Proof of Theorem [3} By the second part of the previous 
lemma, we only need to consider T 2 = 0,l,...,f> — 1. We first 
consider the case when b g is between 1 and p — q — 1. We 
further consider two subcases. 

(i) < r 2 < q: 

Suppose that ( 1221 has no solution for < x < t 2 . As the 
indicator function in (l22l is zero for x = t 2 , t 2 + 1, . . . ,q — l, 
(l22l is reduced to 

x = a g (ri,T 2 ) mod p. 

The number of integers in {t 2 , t 2 + 1, . . . , q — 1}, say d, 
satisfies [d/p\ = m. By Lemma[T| we have either morm+1 
solutions to (f22l for x > r 2 . 

Secondly, suppose that (1221 has exactly one solution for 
< x < r 2 . The indicator function in (f22l is equal to 1 for 
< x < t 2 . Hence, 

0<a g (n,T2)+b g <T2. (27) 

We claim that (l22l has no solution for x — t 2 ,t 2 + 1, . . . , q—1. 
Otherwise, we have 

t 2 < a a (ri,r 2 ) < q, 

which, after combining with the assumption that 1 < b g < 
p — q — 1, yields 

t 2 < a 9 (ri,r 2 ) + b g < p - 1. 

This contradicts with ( fZTI i and proves the claim. For 

q < x < q, 

there are exactly m solutions by Lemma Q] The total number 
of solutions to (l22l for x = 0, 1, . . . , q — 1, is thus rn + 1. 
Hence H 3 i(ti, t 2 ) = m + 1. 

(ii) q<r 2 <p: 

By Lemma [T] (l22l has either or 1 solution for < x < r 2 , 
and either to — 1 or to solutions for r 2 < a; < q. Hence, 
Hgi{i~i,T 2 ) is within the range of {to — 1, m, m + 1}. 

For b g = p — q + 1, . . . ,p — 1, we again consider two 
subcases. 



(i) < r 2 < q: 

By Lemma [T] (f22l has either or 1 solution for < x < r 2 , 
and either to or to + 1 solutions for r 2 < x < q. Therefore, 

H g i(T 1: T 2 ) G {m, m + 1, m + 2}. 

(ii) q<r 2 <p: 

Suppose that (f22l has no solution for < x < r 2 , i.e., 

t 2 < a g (T 1 ,T 2 ) + bg < p. (28) 

We claim that (l22l must have one solution for x in the 
following range 

t 2 < x < p + q. (29) 
From the assumption of q < r 2 < p, we deduce that 

q<p + q-T 2 <p, 

so that the range in ( |29l is non-empty and consists of no 
more than p integers. If the claim were false, we would have 
no solution to (l22l for r 2 < x < p + q, implying that 

q < a g (ri,T 2 ) < t 2 . (30) 

Here, we have used the fact that the indicator function in (l22l 
is equal to zero for x in the range in ( |29l . By adding d30b to 

p-q + l<b g <p-l 

and reducing mod p, we obtain 

1 < a 9 (ri,r 2 ) + b g < t 2 , 

which is a contradiction to ([28). Thus, the claim is proved. 
For x = p + q,p + q + l,...,q — 1, there are exactly to — 1 
solutions to (l22l by Lemma [T] Totally there are m solutions, 
and thus H g i(Ti,T 2 ) = to. 

Finally suppose that (f22l has exactly one solution for < 
x < t 2 . As the number of solutions to (f22l for x = r 2 , t 2 + 
1, . . . , q— 1 is either m — 1 or m by Lemma[T] the total number 
of solutions to d22l is either morm + 1. 

In any case, we see that H g i(ri, r 2 ) is either to, to + 1 or 
to + 2. ■ 

Appendix B 
Proof of Theorem"?] 

We use the following property of Hamming cross- 
correlation which holds for any binary sequence set in gen- 
eral |fT5l 

Lemma 3. Let a(t) and b(t) be binary sequences with period 
L and Hamming weight q. Then 

Y J H ab {r)=q\ 

T = 

We include the short proof for completeness. 
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Proof: 

L-l L-lL-1 

r=0 r=0 t=0 

L-l L-l 

= 5>(*)E & (*- r ) 

t=0 T=0 

L-l L-l 

= J>(t) J>(r) = g 2 . 

t=0 r=0 

The last equality follows from the assumption that the Ham- 
ming weights of a(t) and b(t) are both q. ■ 

1) By Theorem |2j No(j) is nonzero only when j — to or 
j = to + 1, whence, 

N (m)+N (m + 1) = pq. (31) 

On the other hand, we have 

mN (m) + (m + l)N Q (m + l) = g 2 (32) 

from Lemma [3] We can eliminate N (to) from ( BTT i and (|32l 
and obtain 

mpq + No(m + 1) = q 2 , 

which implies jVo(m+l) = q{q — mp) = qq. Substituting this 
into OTb . we get No(m) = pq — No(m + 1) = (p — q)q. This 
proves the first part of Theorem |4] 

2) Suppose that b g = 1, 2, . . . ,p — q — 1. We can set up a 
system of two linear equations in three variables N g (m — 1), 
iV g (m) and N g (m + 1): 

m+l 

E N 3 ( k )=PQ 

k—rn — 1 
m+l 

£ fciV 9 (fc) = g 2 . 

fe=m— 1 

The second equality is due to Lemma [3] Solving for N g (m) 
and N g (m + 1) in terms of N g (m — 1), we obtain © to ||9). 
It remains to evaluate N g {m — 1). The proof is completed by 
showing the following two claims: 

(i) For each k = 0, 1, . . . , to — 1, there are exactly 

b g (p-b g -q) (33) 

order pairs (ri,T2), with < ti < p and T2 = kp,kp + 
1, . . . , (k + l)p — 1, such that H g i{T\,T2) = m — 1. 

(ii) For T2 = top, nip + 1, . . . , g — 1, H g i(ri, T2) does not 
equal m — 1 for all < t\ < p. 

Multiplying d33l by to, we obtain 

N g {m - 1) = m6 g (p - 6 9 - g). 

In the following, we complete the proof of part 2 in the 
theorem by showing (i) and (ii). 

By Lemma [2] part [2] we notice that -ff s i(ri,T2) depends on 
T2 only through the residue of t-i mod p. Hence it is sufficient 
prove claim (i) for k = 0. 

Consider T2 from to p — 1. We want to count the number 
of times that Hgi(ri,T2) = to — 1, for < t\,ti < p. We 
have shown in Lemma[2]that H g i{ri,T2) can be computed by 



counting the number of solutions to (l22l for x between and 
q — 1. Partition the range of x into two disjoint subsets 

Xi := {0, 1, . . . ,p + q - 1}, and 
X2 ■= {p + q,p + q + l,...,q- 1}, 

and consider the number of solutions to (l22l for a; in A'i 
and A2 separately. Because < < p, the indicator function 
1(0 < x < T2) in (l22l is identically equal to for 1 6 A2. 
The number of solutions to (l22l for x G is exactly m — 1 
by Lemma [TJ The problem reduces to counting the number 
of pairs (ti,T2) G Z 2 such that (l22l has no solution for 
x G X\. Observe that a g depends on n and T2 through their 
difference t\ — T2 mod p. We make a change of variables 
( T i: T 2) (u, T2) by defining it as 

7i:=(<7-l)- 1 (r 1 -T 2 )modp. (34) 

Our objective now is to count the number of ordered pairs 
(u, T2) G Z 2 such that the equation 

x = u + b g I(0 < x < T2) mod p (35) 

has no solution for x € X\. We note that b g does not depend 
on T2 and u. 
If 

< u < q, (36) 

then d35l l has at least one solution over x G X\, namely x = 
(u mod p) + p. Indeed, as T2 < p < x, the indicator function 
1(0 < x < T2) is evaluated to 0, and thus if we put x = 
(u mod p) + p in ( 135b , we have (u mod p) + p = u mod p, 
which obviously holds. 
On the other hand, if 

p-b g <u <p, (37) 

then (|35T > also has at least one solution over x G A\ no matter 
what T2 is. Indeed, for < < u, we can set a; = u. Then 
1(0 < x < T2) equals 0, and we see that x = u is a solution 
to (l35l . For u < T2 < p, we can set x = u + b g — p. Then 
we have x < u < and whence 1(0 < x < T2) = 1. When 
u < T2 < p and x = u + b g — p, (l35T l becomes 

u + bg — p = u + b g mod p. 

From (l36l l and (|37] |. H g i(ri, T2) is equal to m — 1 only when 
it = g, g + 1, . . . ,p — b g — 1, For each such u, we now count 
the number of T2 G 1 V for which (l35T l has no solution over 
x G X\. If a; is a solution of d35l l, then a; can take only two 
values, namely u or it + b g . When x = u is a solution, x must 
satisfy x > T2; when x = u + b g is a solution, a; must satisfy 
< x < T2- So, d35l l has no solution over x G X\ if and only 
if for all x G A4, we have u < T2 and u + b g > t%. Putting 
these two inequalities together, we obtain u < T2 < u + b g . 
Consequently, for each u = q, q + 1, . . . ,p — b g — 1, there are 
exactly b g values of T2 such that H 9 i(ri,T2) = m — 1. This 
proves claim (i). 

For claim (ii), we count the solutions to (l35l l for < XT2. 
Since 1(0 < x < T2) is identically equal to 1, by Lemma [Tj 
the number of solutions to (|35T > over < x < T2 is at least to. 
So H g i(ri,T2) cannot be m — 1. 
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This ends the proof of the second part of Theorem |4] 

3) We set up the following system of linear equations in 
variables N g (m), N g (m + 1) and N g (m + 2): 

m+2 

J2N g (k)=pq 

k—m 
m+2 

J2 kN g (k) = q\ 

k=m 

After solving for N g (m) and N g (m+1) in terms of N g (m+2), 
we obtain ( fTOb to ( TlZb . So, we just need to evaluate N g (m+2). 

The evaluation of N g (m + 2) relies on the following two 
claims: 

(i) For each k — 0,1,2, ... ,m, there are exactly 

(p - b g )(q + b g - p) (38) 

ordered pairs (n, T2), with < ri < p and fcp < t 2 < kp + q, 
such that H g i(Ti,T2) = m + 2. 

(ii) For fc = 0, 1, . . . , m— 1, none of the ordered pair (n , r 2 ) 
with < Ti < p and kp + q < t 2 < (k + l)p, satisfies 
H g i{T 1 ,T 2 ) =m + 2. 

Since (i) and (ii) exhaust all possible (ti,T2) € G Pig , we 
multiply (l38l by (m + 1) and obtain 

7V 3 (m + 2) = (m + l)(p-6 3 )(q + 6 3 -p). 

We prove case (i) and (ii) in the rest of this appendix. 

By the second part of Lemma [2] it is sufficient to prove 
case (i) for k = 0. Consider r 2 from to q — 1. Divide the 
range of x into two disjoint parts: 

Af 3 :={0,l,...,g-1}, and 
X 4 := {q,q + l,...,q-l}. 

As in the proof of the second part of this theorem, we make a 
change of variable (t%, T2) — > (u, r 2 ) by defining u as in ( l34l . 
It is noted that A4 consists of mp consecutive integers, and 
the indicator function 1(0 < x < t 2 ) in d35l l equals for all 
x £ A4. By Lemma Q] there are exactly m solutions to (|35T > 
for x £ A4. So, H g i(ri, r 2 ) equals m + 2 if and only if (l35T l 
has exactly two solutions over x £ A3. It reduces the problem 
to counting the number of pairs (u,r 2 ), with < u < p and 
< r 2 < q, such that (f35T > has exactly two solutions over 

X £ A3. 

The only two candidate solutions for (l35l l are x = u and 
x = u + b g . We investigate under what condition both u and 
it + b g are indeed solutions to (l35l l. 21 = it is a solution only 
if 1(0 < x < r 2 ) =0. This implies that it < r 2 . x = u + b g 
is a solution only if 1(0 < a; < t 2 ) = 1. Hence u + b g < r 2 . 
Combining these two conditions, we obtain 

u + bg < r 2 < u. 

This is possible only if u + b g > p, and thus after reduction 
mod p, we have (u + b g mod p) < u. The first necessary 
condition for both u and u + b g are solutions to (l35l l is 

p-b g <u<p. (39) 



Secondly, as it is required that x £ A3, we must have u £ 
A 3 , i.e., 

< u < q (40) 

Putting (|39l and (l40l > together, we have the following 
necessary condition on u. 

p-b g <u<q. (41) 

We note that the range of u in fiTt is nonempty, because 
p — 6 9 < g by assumption, u can take on any of the (b g + q—p) 
value in WH , and for each such u, r 2 can assume values in 
u — {0, 1, . . . ,p — b g — 1}. Hence, the total number of pairs 
( T i: r 2) such that < r 2 < g and -ff 9 i(ri,r 2 ) = m + 2 is 

(p-bg)(bg - (p-q)). 

For case (ii), we again use the fact that H g i(ri, r 2 ) depends 
on r 2 only through the residue of r 2 mod p, and establish 
case (ii) only for k — 0. Let r 2 be in the range q < r 2 < p. 
Consider the solutions to d35l l separately in < x < r 2 and 
r 2 < x < q. For < x < r 2 , 1(0 < x < r 2 ) is identically 
equal to 1 . By Lemma Q] there are at most 2 solutions to (|35| > 
for a; = 0, 1, . . . , t 2 — 1. For t 2 < x < g, 1(0 < x < r 2 ) 
is identically equal to 0, and by Lemma Q] there are at most 
m — 1 solutions to d35l l for r 2 < x < q. There are totally at 
most m + 1 solutions to (|35| >. 

This completes the proof of Theorem |4] 

Appendix C 
Proof of Theorem[9] 

In this proof, p is a prime number and q is an integer 
relatively prime to p and strictly larger than 2p 2 , L = pq 
is the sequence period, and 7 is the multiplicative inverse of 
p mod q, i.e., 7p = 1 mod q. The unique integer between 
and q — 1 which is equal to x mod g is denoted by (2 mod q). 
The translate of a subset S in G p ^ q by (ti,t 2 ) is defined as 

S + (t 1 ,t 2 ) :={{x,y) + {T 1 ,T 2 ) : (x,y) eS}. 

Consider user g, where i = l,2,...,p — 1. If user g starts 
transmitting at time to, then the channel-activity signal is 
matched to s g (t) at time to. The receiver will never fail to 
detect the presence of user g, meaning that if user g does 
start transmitting, the receiver can always detect this change 
of status from idle to active. The only sources of error are 
(a) detecting a user but in fact that user is not active, and (b) 
miscalculation of the start time. We refer to the error in (a) as 
false alarm and (b) as synchronization error. 

We now show that false alarm cannot occur. Suppose on the 
contrary that the channel-activity signal is matched to s g (t) at 
time to, but user g is idle from time to to to + L — 1, If this 
happened, the q time indices in I g + to would be covered 
by the protocol sequences of other active users. However, 
the cross-correlation between user g and each other active 
user is upper bounded by [q/p\ + 2, by Theorem [3] Because 
user g is assumed to be inactive in this period, the number 
of simultaneously active users does not exceed the maximum 
(p + l)/2, and hence the number of time slots in I g + to with 
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a packet or collision observed is no larger than 

'p+l 



(L<z/pJ 



/l 2xp+l 
= q \p + - q ) — 
fl 2 \p+l 

<q {p + 2y) — 

_ q ip + 1\ 2 
~ 2\~~p~ 



(42) 



In (l42l . we have used the assumption that q > 2p 2 . We note 
that for all p > 3, the factor (p + l) 2 /p 2 is strictly less than 
two. We obtain 



(L^/pJ 



But g is precisely the total number of ones in a period of s g (t). 
The q time slots indices by I g + t + cannot be covered by 
any other (p + l)/2 CRT protocol sequences. The channel- 
activity signal c(t) cannot be matched to s g (t), and therefore 
false alarm cannot occur. 

For synchronization error, assume that user g is idle from 
time to — L + 1 to to — 1, and becomes active at time io- 
Our objective is to show that the channel-activity signal is not 
matched to s g (t) at time to — t, for any integer r between 
1 and L — 1. The idea of showing that synchronization error 
cannot occur is the following. If the channel-activity signal 
were matched incorrectly to s g (t) at to — r, then the receiver 
would observe q "1" or "*" at time slots indexed by to — 
t + ^'~Hl g:Ptq ). Among these q time slots, say b of them 
come from a shifted version of s g (t), starting at time to. We 
then show that the remaining q — b slots cannot be covered 
by the other active users. We divide the proof into several 
propositions below. 

In the modified CRT correspondence an element t € Z pq 
is mapped to (t modp, yt mod q). In order to visualize the 
mapping, we introduce a matrix M p 9 . 

Definition 8. Given relatively prime integers p and q, let 
M p>g be a p x q matrix whose (i,j)-entry equals t if 
i = t mod p and j = jt mod q, for t = 0, 1, . . . ,pq — 1. The 
rows and columns of M p g are indexed by {0, 1, . . . ,p — 1} 
and {0, 1, . . . , q — 1}, respectively. 

Each integer from to pq — 1 appears exactly once in M Pi q. 
An example for p = 3 and q = 19 is shown in Fig. |4] We pay 
special attention to the integers from to p 2 — 1, and want to 
get a handle on where they are located in M p g . In Fig. [4] we 
can see that 0, 3 and 6 are on the upper left corner of M 3 lg . 
The numbers 1, 4 and 7 occupy three consecutive entries in 
row 1. The numbers 2, 5 and 8 occupy three consecutive 
entries in row 2. 



Proposition 1. ( i) For any i and j, the l)-entry is equal 

to p plus the (i, j)-entry mod pq. 

( ii) Under the modified CRT correspondence $ p , the 
integers kp, for k = 0, 1, 2, . . . , q — 1, are mapped to (0, k). 
They appear in the first row q/M Pj g. 

( Hi) The numbers from 1 to p 2 — 1, except the multiples of 
p, are located between column 2p + 1 and column q — p — 1 
inclusively in M p ( j. 

Proof: (i) The (0, l)-entry in M p q is labeled by p, 
because 

p h- > (p mod p, 7J? mod q) = (0, 1), 

and 7p = 1 mod q by the defining property of 7. 

(ii) For k = 0,1,2,..., q- 1, 

^ q {kp) = (kp,k P1 ) = (0,fc(l)) = (0,fc). 

(iii) We have the following claim: 



2p < (fc7 mod q) < q — 2p 



(43) 



for k = 1, 2, . . . ,p — 1. 

We prove the claim by contradiction. Suppose that kj, after 
reduction mod q, is between 1 and 2p. Then, k'jp mod q is 
equal to p, 2p, 3p, . . . , or 2p 2 . Since q > 2p 2 , these p numbers 
remain unchanged after reduction mod q. However, k^p = 
k(^p) = k mod q, and this contradicts the assumption that k 
is between 1 and p— 1. Now suppose that (^7 mod q) is equal 
to q — 2p + 1, q — 2p + 2, . . . , or q — 1. Then, the value of 
kjp, after reduction mod q, is equal to 

q - 2p 2 + p, q - 2p 2 + 2p, . . . , or q - p. 

Since k^p = k mod q, this also contradicts that k is between 
1 and p — 1. This finishes the proof of the claim. 

Let £ be an integer between 1 and p 2 — 1 which is not a 
multiple of p. We can write £ as mp + k for some m and fc 
between 1 and p — 1. By part (i), the location of ^ in M p 9 
is m steps to the right of the location of k in M p r But fc 
cannot be located to the right of column q — 2p. The right-most 
column in M p g which may contain £ is thus q — p — 1. This 
finishes the proof of part (ii). ■ 

Proposition 2. (/) Let S = {0, 1, 2, . . . ,p 2 — 1}, and S' be 

the image of S under $ p q , i.e., S' = $ p q {S). We have 



\l g r\(l h + (n,T 2 ))nS'\ <2 



(44) 



for any given (n, T2) ant/ (? 7^ /i. 

(7j) For g = 1, 2, . . . ,p — 1, f/zere are exactly p ones in 
the first p 2 bits of s g (t), i.e., there are exactly p ones among 
s 9 (0), Sff (l),..., Sff (p 2 -l). 
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The quantity on the left hand side of ( |44T > can be interpreted 
as the partial Hamming cross-correlation, defined as 

p 2 -i 

J2 s g (t)s h (t + r). 

t=Q 

We only consider the number of overlaps in the first p 2 
time indices. The proposition asserts that the partial Hamming 
cross-correlation of two CRT sequences cannot exceed two. 

Proof: (i) By part (i) of Prop. Q] for each k between 
and p — 1, the following p integers, 

k,k + p,k + 2p, . . . ,k + (p - l)p, 

occupy p consecutive horizontal entries in the fcth row of 
M p g . Wrapping around the right boundary of M p g is pre- 
cluded by Prop. [T] 

A common element of X g and Xh + (ti, t 2 ) is in the form 

(gt 1 ,t 1 ) = (ht 2 +n,t 2 + T 2 ) (45) 

for some t\ and t 2 between and q — 1. If t 2 + t 2 is between 
and q — 1, then the first coordinates of the two order pairs 
in (|45t are equal to ht 2 + t 2 , with t 2 satisfying 

g(h + t 2 ) = ht 2 + ti mod p. (46) 

If t 2 + t 2 is larger than q, then t% = t 2 + t 2 — q, and the first 
coordinates of the two order pairs in ( l45T > are equal to h 2 + t 2 , 
with t 2 satisfying 

g(t 2 +t 2 - q) = ht 2 + Ti mod p. (47) 

Hence, ht 2 + t\ may assume only two values mod p, one 
from ( |46*| > and the other one from (1471 . Let A g be the p x q 
array with characteristic set I g . We see that the elements in 
I g n (Z/j + (ti,t 2 )) are located in at most two rows in A g . 
Since p consecutive entries in a row of A g contain exactly one 
"1", at most two elements in I g n {Ih + (tx,t 2 )) are covered 
by S'. 

(ii) Let A g be the p x q array with X g as the characteristic 
set. From the remark before Definition |4] any p consecu- 
tive columns in A g form a permutation matrix. For each 
b = 0, 1, . . . ,p — 1, the time indices 

b, b+p,...,b+(p-l)p 

are p consecutive entries in a row in M p 9 . Hence there is 
exactly one "1" among s g (b), s g (b + p), . . . , s g (b + (p— l)p). 
Since this is true for b = 0, 1, ... ,p — 1, we conclude that 
there are exactly p "1" in s g (0), s g (l), . . . , s g {p 2 — 1). ■ 

Proposition 3. For y = 0, 1, . . . , q — p, let V' be the index 
set 

V' y :={(i,i)eG M : o<i<p-i,y<j <y+p-i}- 

Let V y be the corresponding set of time indices in Z pg under 
the modified CRT correspondence, 

Then 

(i) For g 7^ h, and delay offset (n, t 2 ), 

\i g n(i h + (T 1 , T2 ))nr' y \<2. 



(ii) For g — 1,2, ... ,p, there are exactly p ones in s g (t) 
for t e V v . 

Proof: The elements in V' can be regarded as a square 
submatrix in a p x q matrix. The proof of part (i) of Prop. [3] 
is similar to that of Prop. [2] and is omitted. The second part 
follows from the fact that every p consecutive columns in A g , 
which stands for the p x q matrix with characteristic set I g , 
form a permutation matrix, and hence contains exactly p ones. 

■ 

Proposition 4. Let S', V' y be defined as in Prop.\2\and\3] For 

g = 1, 2, . . . , p — 1, and r between 1 and L — 1, the subset 
B g ^ of time indices t in {0, 1, 2, . . . , L — 1}, which satisfy 

Sg (t) = l/or* = 0,l,...,T-l, 
s g (t ) — 1 and s g (t — t) = for t = r, r + 1, . . . , L — 1, 

contains ^q{S' nig), or &p,q{V y C\l g ) for some y. 

Proof: Let s g (t) be the acyclic shift of s g (t) to the right 
by delay offset r. The time indices in B 9tT correspond to the 
"1" in s g (t) which is not covered by s g (t). We consider two 
cases: (i) p 2 < t < L, and (ii) 1 < r < p 2 . 

Case (i). When p 2 < t < L, the first p 2 bits in s g (t) are not 
covered by s g (t), and therefore B g T contains ^'~^(S' C]I g ). 

Case (ii). Recall that in the proof of Theorem [6] it is 
mentioned that the intersection of I g and I g + (ti ,t 2 ) is either 
empty, or an arithmetic progressions in G Pi9 with common 
difference (5,1). The intersection is non-empty if and only if 
(t~i,t 2 ) equals (g, l)k for some k = ±1, ±2, . . . , ±(q — 1). 

Suppose that r is a nonzero multiple of p between 1 and 
p 2 - 1. By part (i) of Prop.Q] &p iq (r) = (0, r/p). (r is one of 
the p left-most entries in the first row of A4 p , q .) In this case, 
$'(t) does not equal (g, l)k for any k 6 {1,2, . . . ,p — 1}. 
(We have used the assumption that g is non-zero in this step.) 
This implies that the ones in A g with indices in V' T , p are not 
covered by s g (t). Thus, 

Now suppose that r is between 1 and p 2 — 1 but is not a 
multiple of p. Let (r l5 r 2 ) = $' p q (r). Using similar argument 
as in the previous paragraph, if (ti, t 2 ) tfL1 g , then 

6 9 , T D$^(?; 2 ni s ). 

Otherwise if (ti,t 2 ) ^ I g , then the intersection of X g and 
X g + (n,r 2 ) equals 

{(g,l)k: fc = r 2 ,r 2 + l,...,g-l}. 

Consider the p columns in M p g to the left of r, namely the 
time indices associated with V T2 _ p in M p g . The correspond- 
ing time slots are not covered by s g (t). Therefore, 

■ 

Suppose that s g (t) is actually transmitted at time t and 
the receiver tries to match the channel-activity signal c(t) 
with s g (t) at time t — t (See Fig. [5]). If the channel-activity 
signal were mistakenly matched to s g (t) at to — t, then by the 
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Fig. 5. The protocol sequence s g (t) is transmitted at time 4q. Sequences 
si(t) and S2(t) are transmitted at some other time. The receiver tries to match 
c(t) with s g (t) at to — T. 



previous proposition, either (i) the p "1" in the first p 2 bits of 
s g (t), which occur at time indices 

($^ 1 (i 3 )n{o,i,...,P 2 -i}) + t-r, 

or (ii) the p "1" in s g (t), which occur at time indices in 

for some y between and q — p, are covered by the other 
active users. By Prop. [2] and (O, each of the other active users 
can contribute at most two overlapping slots. As there are no 
more than (p — l)/2 other active users, the total number of 
ones that can be covered by other active users are at most 
2 • (p — l)/2, which is strictly smaller than p. This proves 
that the channel-activity signal cannot be matched to s g (t) at 
to — t for any r = 1, 2, . . . , L — 1. This completes the proof 
of Theorem [9] 
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